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Recently it was found that the dynamics in a Heisenberg spin-chain subjected to a 
sequence of periodic pulses from an external, parabolic, magnetic field can have a close 
correspondence with the quantum kicked rotor (QKR). The QKR is a key paradigm of 
O ■ quantum chaos; it has as its classical limit the well-known Standard Map. It was found that 

a single spin excitation could be converted into a pair of non- dispersive, counter-propagating 
spin coherent states equivalent to the accelerator modes of the Standard Map. Here we 
consider how other types of quantum chaotic systems such as a double-kicked quantum rotor 
or a quantum rotor with a double-well potential might be realized with spin chains; we 
discuss the possibilities regarding manipulation of the one-magnon spin waves. 

§1. Introduction 

It has been shown that the quantum properties of certain many body systems 
may be analyzed by considering the dynamics of an analogous 'image' one-body 
system. 1 -* In the case of certain types of spin-chains, there can be one-body image 
systems with a well-defined classical limit, which can be chaotic or integrable. 2 ) 
*^ In another context, there has also been considerable interest in the dynamics 

of quantum spin-chains because of their potential for quantum information appli- 
cations. Quantum state transfer is one example: the ability to transfer a qubit, or 
some arbitrary quantum state, with high fidelity along a spin-chain has been ad- 
dressed in several works; in 3 ) quantum state transmission of a single spin-flip along 
a ferromagnetic Heisenberg chain was investigated. In a subsequent work 4 ) it was 
shown that such a chain, in the presence of an external, static, parabolic magnetic 
field, can give perfect transmission if the initial spin state is a specific coherent state. 
Recently 5 ) a close correspondence between the unitary time evolution of the 



ferromagnetic Heisenberg spin-chain -for a single spin-excitation- and that of the 
Quantum Kicked Rotor (QKR) was noted. To complete this correspondence, a 
periodically-pulsed parabolic magnetic field must also be applied. The QKR is a lead- 
ing paradigm of quantum chaos 6 ) which has been extensively studied theoretically. 
Its classical counterpart is one of the best-known 'textbook' examples of Hamiltonian 
chaos, the Standard Map. 7 ) The QKR has also been investigated experimentally- 
since it was shown by a cold atom group in Texas that a realization with cold atoms 
with pulsed standing waves of light is possible. 8 ) A number of other aspects of this 
system in its quantum chaotic regime were also studied by this and several other 
experimental cold atom groups. 9 )' 10 ) 

In 5 ) it was found that with the pulsed parabolic field, one can employ certain 
expressions found for the QKR and the Standard Map to describe the entanglement 
properties. It means also that the system exhibits not only generic forms of quan- 
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turn chaotic behaviour in the spin-chain like exponential localization (analogous to 
Anderson localization seen in disordered metals) but also generate phenomena such 
as 'accelerator modes' which are a feature specific to the Standard Map alone. 

In £j2] we review the QKR one-body image of the ferromagnetic Heisenberg spin- 
chain. In §3 we show that the double-kicked rotor (QKR-2), a system which has 
been the subject of an experimental study 10 ) and has been found to have rather 
different dynamics from the QKR, also has a spin-chain analogue. We propose that 
rather than implementing the QKR-2 which is a KAM dynamical system, it may be 
easier to implement a closely related random map. In §4 we show that double- well 
kicked rotors may be implemented with spin-ladders or ferromagnets with next- 
to-nearest neighbor exchange interaction. To conclude, we discuss the possibility 
of investigating the corresponding antiferromagnetic dynamics and we discuss the 
potential uses of the quantum chaos as a means to manipulate spin waves. 

§2. Quantum kicked rotors with Heisenberg ferromagnetic spin-chains 

The Heisenberg spin-chain has the well-studied Hamiltonian: 

H hc = -J^ n ^ n+l -Y. B ^- W 

n n 

For a ferromagnet, J > 0. We consider the effect of applying a periodic sequence 
of short pulses from a parabolic magnetic field. We can model these by a series of 
8— kicks. Note that in equivalent QKR atomic experiments the pulses will, of course, 
have finite duration. The combined Hamiltonian is: 

N B 

H = Hhc + ^2 ~2~ ( n ~~ no) 2 "? £ *(* " J" r o) (2) 

n=l j 

where To is the period of the pulses; Bq is the amplitude of the applied parabolic 
magnetic field; the length of the chain, N > 100 in the present work. In 5 ) this 
effect was investigated for a single excitation on an open chain. Here we consider the 
simpler case of a chain with periodic boundary conditions (eg. a ring). However, for 
long chains, the dynamics we explore are not sensitive to these boundary conditions. 

Since [H, S z ] = 0, where S z is the total spin component in the z direction, if 
we prepare an initial quantum state comprising all spins pointing 'up' and a single 
spin 'down' at some arbitrary site ip(t = 0) = |so), one may, for all t, obtain ip(t) in 
terms of a basis of states |s), which have a spin-down at a single site s on the chain 
but all other spins up (along z). The eigenstates of Hh c in this basis distribute the 
spin-flip periodically along the chain: 

1 N 

|m) = -= (3) 
V j = i 

These states represent spin waves (magnons) with wavenumber k m = 2m7r/N, 
2m 6 (-N, N] and may be obtained using the Bethe ansatz. The corresponding 
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eigenenergies are given by the one-magnon dispersion relation for a ferromagnet: 

E m -E = J(l - cos km) + B, (4) 

where, Eq is the ground state energy —JN/4 — NB/2. 

An analytical form for the time-evolution operator f//j C (t, 0) = exp{— iHhct] 
may be given, in the single-excitation basis, as a matrix of elements: 



U hc (To) = e~ i ~&^ r ~ n °) 2 J_ e i ( r - s ) k m+iJT cos k m ^ 

Here, we disregard the overall phase due to Eq and the uniform static field B (or 
formally set 2BTq = 2ir). In an actual realization, the external field is significant 
as it introduces a gap with the fully aligned ground state: it does not modify the 
dynamics of interest here. 

In comparison, the QKR has Hamiltonian: 

p2 

H = — - K cos x^5{t-nT). (6) 

n 

Free evolution of the particles is followed by a short kick from a sinusoidal potential. 
K (stochasticity parameter of the Standard Map) is the kick-strength; it is related 
to the intensity of the optical lattice in atomic experiments. In a basis of momentum 
states \l), where p = lh, the QKR unitary evolution operator has matrix elements: 

U% KR (T) = e~ i2 T 12 — ^ e *(n-l)x+if cos x ^ (7) 

From ([5]) it is evident that the 'free evolution' part of U due the static spin- 
exchange interaction is simply a discretized version of the kick part of JJ^ KR ; while, 
in turn, the 'kick' provided by the parabolic magnetic field is equivalent to the free- 
evolution of the QKR. The magnon wavenumber maps onto position for the QKR, 
ie k m — ► x; similarly, the spin site maps onto momentum, so s — > I. Since, for the 
QKR, we rescale time units so T = 1, we see that in ([5]), Bq — > H ie the parabolic 
field can be considered an effective value of h in the QKR image system. Finally, 
the classical stochasticity parameter, K = JTqBq. 

Equations © and ([7]) are sufficiently close that the behavior of the QKR maps 
onto the Heisenberg spin chain at even the detailed level. We see, not only, expo- 
nential localization for strong K, but also, for values of K slightly above K = 2ir, 
we see accelerator modes. Accelerator modes are a peculiarity of the standard map: 
deep into the full chaotic regime, small transporting stable islands re-appear and 
cause ballistic rather than diffusive transport over part of phase-space. This is il- 
lustrated in Fig[TJ We show, for K ~ 6.6, the effect of repeated application of the 
unitary matrix of © on a state initialized with a spin flipped 'up' on a site for which 
s ~ 2ttM/Bq where M = 0, 1, 2, ... which is close to an accelerator island. In Fig{T] 
we took simply M = and ip(t = 0) = |no); ie a site near the center of the chain. 
The spin-amplitude spreads out into a very irregular 'chaotic' distribution around 
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Fig. 1. Effect of accelerator modes in a Heisenberg spin-chain. P(s) represents the probability of 
finding the excitation at site s. The accelerator modes are the 'spikes' at the leading edge of 
the distribution. They correspond to a counter-propagating pair of coherent states ejected from 
the centre. We take Bq = 1/15 and JT = 100. The upper line is at t = 3Tb; the lower curves 
correspond to consecutive periods jTo with j = 4, 5, 6... as numbered. The dotted line indicates 
the form of the parabolic field (scaled by a constant factor) which is pulsed on/off every period 
at times t = jTo. The accelerator modes represent over 25% of the total probability; they 
advance an equal distance (shown below to be 2-k/Bq ~ 94 spin sites) each period, and after 
just 3 pulses are well separated from the central, 'chaotic' remnant. 



the site n$. But most strikingly, we can see a pair of counter-propagating spikes, 
'hopping' around 2tt/Bq ~ 94 spin-sites each consecutive period. 

Away from parameters where there are accelerator modes, but still in the chaotic 
regime K > 4, the QKR exhibits an analogue of Anderson localization in disordered 
metals; for the QKR, this means exponential localization of momentum. For the 
Heisenberg spin-chain one can straightforwardly infer (and verify numerically 5 )) that 
a single excitation will, after a timescale t > t* ~ (JTo) 2 give rise to an exponential 
spin probability distribution: 

P( S )~|exp{-2| S - So |/^}- (8) 

where L ~ (JTo) 2 /4. Note that the variance of the spin distribution < (s — sq) 2 >— 
L 2 is independent of the initial state \sq). Similarly, the asymptotic variance of the 
QKR momentum distribution < (p — Pq) 2 >~ K 2 /h is independent of the initial 
momentum of the atoms, since at t = 0, typically, the momenta of the cold atom 
cloud is close to some value Pq. Typically N(p) ~ exp (p — Pq) 2 /a 2 , where a ~ 4 in 
the usual units depends on the initial temperature of the cold atoms . 
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§3. Double-kicked rotors (QKR-2) with spin-chains 

In 10 ^ the double-kicked rotor (a variant of the QKR) was investigated experi- 
mentally with cold atoms: a cloud of Cesium atoms was subjected to pairs of closely 
spaced pulses from standing waves of light. The corresponding Hamiltonian is: 



p2 

H = — K cos x 



^5(t - nT) + 5{t - nT + e) 



0) 



with period T, kick strength K and e <C T the time between each kick in the pair. 

The variances of the asymptotic (long time) distribution in the QKR-2 depend 
strongly on the initial state as shown in Figj2j The results display 'momentum 
trapping' regions where atoms absorb little energy if Pq ~ (2m + l)vr/e. Here the 
variances remain close to the initial values (which depend on the temperature -of 
order fiK- of the cloud. In between these trapping regions, the typical variances 
~ 7r/e are, unlike the QKR, independent of both the kick strength and of h. 

Fig|2] also shows a typical surface of section. Classical phase-space shows a 
cellular structure, with the cell-boundaries corresponding to the momentum trapping 
regions. In 11 ) properties of the QKR-2 were investigated. For T ~ 1 — 2, the 
trapping regions correspond to thin mixed-phase-space regions filled with islands and 
cantori. As T — > oo all islands disappear and a non-KAM limit, with no islands, is 
approached. But a remarkable feature of this system's dynamics is that the trapping 
effects remain as strong, even as T — > oo. 



300 




Fig. 2. Left: Experimental momentum variances for the 2-QKR, obtained with cold atoms. ' 
The figure shows that at trapping sites Po ~ (2m + l)7r/e there is little spreading, while in 
between, < (p — Po) 2 >~ {jv/e) 2 . Right: Surface of section plot for the 2<5-KP, K e — 0.35, for 
arbitrary, large r e —* oo, illustrating the classical dynamics of the 'one-body image' of the double 
pulsed spin chain . The cellular structure is evident: phase space is divided into regions of fast 
diffusion separated by porous boundaries, ie narrow trapping regions. For the equivalent spin- 
chain, this would correspond to spin sites s for which sB t ~ ±(2m + l)-7r. where m = 0, 1, 2, ... 
Spin excitations prepared at these sites would remain highly localized, while those prepared in 
between trapped sites would spread over one whole cell. 
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The classical map for the 25-KP is a simple extension of the Standard Map: 

PN+1 = PN ~ K sin xn; Pn+2 = PN+l - K sin xn+i 
xn+i = x N + Piv+ie; x N+2 = x N+ \ + pn+2T 

(10) 

where r = T — e is the time interval between the kick-pairs. Starting from the map 
(|10p with N = we re-scale all variables: p e = pe, K e = Ke and t £ = t/c > 1 to 
obtain a re-scaled map for which the momentum 'cells' are of width 2ir. While it is 
usual for the kicked rotor to be re-scaled in terms of T, the appropriate scaling here 
is in terms of e: 

Pi =Pl~ K e sin x o] P\= Pi ~ K e sin xi 

Xt = XQ+p\] X 2 = Xi+p%T e . (11) 

In the atomic experimental and theoretical studies, 10 )' n ) the regime r e = 20 — 
100 was investigated. However, it was found in 12 ) that diffusion correlations which 
determine the trapping depend only on K e ; the dynamics is completely insensitive 
to t of r e > 50. In effect, very similar dynamics is obtained if we replace xq in 
([lip by a random number in the interval [0, 2tt]. The classical diffusion/trapping 
depends only on correlations between kicks in each pair; consecutive kick pairs are 
completely uncorrelated for r e > 50. We can exploit this useful feature to facilitate 
an implementation of the QKR-2 with a Heisenberg spin-chain. 

The time evolution operator for this system may be written: 

Uqkr-2 = eH^L^f c H e H^] e Hr co H. (12) 

It is easy to see that U is invariant if the products K e = Ke and h e = he are 

kept constant; while the free propagator Jjj Tee = e~ % ^~ simply contributes a near- 
random phase. Provided that l 2 Th 3> 2tt, the results are quite insensitive to the 
magnitude of (T — e)h. 

Now, for the image Heisenberg spin chain of the QKR-2, the quantum map 
corresponding to ([T2"j) is: 

Ul c = e HEL ^(n-n ) 2 *?} e [-iH hc (T )] 

e H£n=l ^(n-n f*?] e [-iH hc (T )]^ ^3) 

Hence evolution under the Heisenberg Hamiltonian H^c for a period To is followed 
by a short impulse from a parabolic magnetic field; however now the strength of the 
parabolic field alternates in value between a strong field B T and a weak field B e . We 
see that the ratio between the fields corresponds to the parameter r e of the QKR-2, 
so B T /B e » 1 = r e . 

If we alternate the field strengths, we find that the image spin-chain will show all 
the dynamics of the QKR-2: a cellular structure for the phase space, with most initial 
states evolving to a final state with variance of the spin-probability < (s — sq) 2 > 
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Fig. 3. Double-kicked spin chain with the strong parabolic field replaced by a randomly varying 
field in time. This effectively removes correlations between 'kick pairs' in the evolution and 
causes spin excitations to remain trapped within one segment of the chain, or between segments 
if prepared in a trapping region. The figure shows two adjacent cells for B t = 0.025 and JTo = 7, 
leading to effective K t — 0.175; so denotes the starting position of the excitation along the chain. 
Excitations that start at the center of the chain (so = 500) spread over one whole cell but no 
further. An excitation starting off-center in the neighboring cell at so = 700, spreads over that 
cell, but not beyond the trapping region at so — 600. Excitations prepared near the center of the 
trapping region remain highly localized over about 30-40 spin sites. Near the edges excitations 
can 'escape', but the probability is very low as shown. 



(tt/ B e ) 2 . In other words, for JTq 3> 1 the single excitation spreads out more or less 
uniformly over a segment of the chain of width ~ ir/B e but remains trapped within 
this segment. States prepared at the edges of the cells would remain highly localized 
over just a few spin sites. 

However, as it is reasonable to suppose that a weak parabolic magnetic field 
might pose less of a technical challenge than a strong parabolic field, we note that 
the strong field can be simply replaced by a random phase, ie in (|13j) we replace the 

term exp —i Y^ n =l %~( n ~ n o) 2(7 z = ex P ~i(PnO~z) where f3 n is a magnetic field of 
randomly varying magnitude, along the z direction. In Figj3]we show spin probability 
distributions obtained for the double spin chain where we take the random field (5 n 
to be uniformly distributed in the interval [0, 2ir]. This corresponds to the dynamics 
of the (non-KAM) r e — > oo limit. 

§4. Double-well kicked rotors and anti-ferromagnets 

Spin-chains in general can have longer ranged exchanged interactions beyond 
the nearest-neighbor type of ([T|). In particular, next-to-nearest neighbor (NNN) 
interactions are well-known and are a feature of spin ladders. 13 ^ 

H hc = E (- J i° n ■ aH+1 - J ^ n ■ an+2 ) - E B < ( 14 ) 



We restrict ourselves now to the case where the leading term is ferromagnetic 
(J\ > 0) but allow the NNN term to have either a ferromagnetic {Ji > 0) or anti- 
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ferromagnetic ( J2 < 0) form. While the latter can lead to magnetic frustration and 
become gapless, we do not concern ourselves with that parameter range. 

A simple Bethe-ansatz treatment of (|14p will obtain similar eigenstates to © 
but a modified dispersion relation: 

E m — Eq = (Ji + J 2 — J\ cos k m — J2 cos 2k m ) + B (15) 

where Eq = — (J% + J2)N/4 — BN/2 for the spin ladder. Hence this system maps 
onto a kicked rotor with a 'double-well' potential: 

pi 

H = — - (Kt cos x + K 2 cos 2x) ^ S(t - nT), (16) 

~ n 

where take the same signs as J\%. This type of potential has slightly different 
dynamics from the QKR; if the kicking period (or for the spin-chain, the quadratic 
field) is varied, it can even give rise to a type of Hamiltonian ratchet. 14 ) 

However, as a potentially useful means of manipulating spin-waves, it is worth 
considering the regimes with low kick strength (small K\,K-i) where the classical 
dynamics has large stable islands. We recall that, in the absence of the parabolic field, 
the Heisenberg ferromagnet (one-magnon) dynamics is equivalent to the resonant 
QKR: quantum states are all fully delocalized and a single spin flip will spread over 
the whole spin-chain. However, kicking by a weak field will generate large stable 
islands where spin-waves are trapped. The location of these islands will determine 
whether spin excitations give rise to low or high energy magnons - low or high k m . 

For simplicity, we consider only the case J\ = J 2 > as well as J\ = — J 2 > 0; 
we note, however, that it may be possible to adjust the relative magnitude of J\ and 
J2 (possibly by varying the distance between the adjacent spin-chains). FigJH shows 
the classical surfaces of section for the ferromagnetic NNN and antiferromagnetic 
NNN respectively. 

Finally, we consider briefly the case of an antiferromagnetic nearest-neighbor 
interaction. The Bethe ansatz treatment obtains a one-magnon dispersion relation 
for the spin waves E m = J\ smk m \ , hence for k ~ the observed antiferromagnetic 
dispersion relation takes a linear E m ~ J|fc m | form. This would make the equivalent 
kicked rotor dynamics somewhat analogous to classical maps with discontinuities 
such as the tent-map or saw-tooth map. Unfortunately, for usual antiferromagnets, 
these states are gapless. The E m = J\ sin k m \ relation represents a lower bound to a 
continuum of excitations. It may be that an NNN type interaction could introduce 
a 'gap' between the linear dispersion relation and the continuum. Nevertheless, the 
prospect of a practical realization of a type of 'saw-tooth map' using spin chains 
would represent a much harder technical challenge. 

§5. Conclusion 

If parabolic fields can be generated one can show that the magnetic field strengths 
required are not unreasonable. For example, we assume the parabolic field ranges 
from B = — 10~ 6 across a spin chain of 10 sites (B = 1CP 6 atomic units cor- 
responds to 0.47 Tesla, a modest laboratory field). Hence Bq ~ 10~ 14 . Since we 
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Fig. 4. Surface of section plots for the double-well kicked rotor with K\ = 0.35 and K2 = ±Ki. 
Stable islands can be used to influence the time evolution of spin waves in a spin ladder with a 
ferromagnetic NN interaction, (a) A ferromagnetic NNN interaction corresponding to K3 ~ K\ 
will mean that a spin-wave packet excited at low fc m will remain confined at low energies; while 
in (b), an antiferromagnetic NNN interaction (K2 = —Ki) will favor intermediate k m . 



require N 2 BQ&t S> 1, where St is the duration of the short magnetic pulse, we 
need St ~ 10 8 — 10 9 au ie St ~ 1 — 10 nanoseconds. We note that picosecond 
pulses are possible in current spin- wave experiments. We also require 2JTq S> 1, 
but for the 'split-operator' approximation in ([5]), we implicitly assume 2 J St <C 1. 
Since J ~ 1GHz for ordinary Heisenberg ferromagnets, this not unreasonable, if 
To ~ lfisec, though it would be better to have a smaller value of J ~ 0.1GHz. The 
smaller J is, the smaller Bq can be (and the longer the duration of the pulses). For 
the double-kicked rotor, provided it is possible to introduce a random phase, the 
parabolic field could be 100 — 1000 times smaller. 

We conclude that if parabolic magnetic fields can be applied to quite common 
ferromagnetic spin chains, the classical chaotic dynamics of the QKR would suggest 
a number of new ways of manipulating the transport of the spin waves. 
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